%% Load data
clear;
close all;
load('OD_data_array_Feb18.mat');

%% Plot initial 25C shifts

ms = 8;

t = A(:,1); % 25C to 37C shift
od = A(:,2);
dod = diff(log(od))./diff(t)*60

figure;
plot(t,log(od), 'ro-', 'markersize', ms, 'MarkerFaceColor','r');
ylim([-2.5 0]);



t = B(:,1); % 25C steady-state
od = B(:,2);

hold on;
plot(t,log(od), 'bo-', 'markersize', ms, 'MarkerFaceColor','b');
ylim([-2.5 0]);
xlim([-30 100]);
set(gca, 'fontsize', 20);
set(gcf, 'position', [0 0 400 300])
xlabel('Time (min)')
ylabel('Log(OD)')
box off;

legend '25 to 37' '25'

dod = diff(log(od))./diff(t)*60


%% Plot initial 37C shifts
t = C(:,1); % 37C steady-state
od = C(:,2);
dod = diff(log(od))./diff(t)*60

figure;
plot(t,log(od), 'ro-', 'markersize', ms, 'MarkerFaceColor','r');
ylim([-2.5 0]);


t = D(:,1); % 37C to 25C shift
od = D(:,2);
dod = diff(log(od))./diff(t)*60

hold on;
plot(t,log(od), 'bo-', 'markersize', ms, 'MarkerFaceColor','b');
ylim([-2.5 0]);

set(gca, 'fontsize', 20);
set(gcf, 'position', [0 0 400 300])
xlabel('Time (min)')
ylabel('Log(OD)')
box off;

ylim([-2.5 0]);
xlim([-30 100]);

legend '37' '37 to 25'

